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. Abstract 

We present a computer code that analytically evaluates the matrix elements of the mi- 



croscopic nuclear Hamiltonian and unity operator between Slater determinants of displaced 
gaussian single particle orbits. Such matrix elements appear in the generator coordinate 



X 

model and the resonating group model versions of the microscopic multicluster calculations 
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PROGRAM SUMMARY 8 P t 

Title of the program: MCKER 



Catalog number: 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland (see applica- 
tion form in this issue) 

Licensing provisions: None 

Computer: IBM PC 386DX 

Operating system: DOS 3.2 

Programming language used: Mathematica v. 1.2, Fortran 

No. of lines in distributed program,including test data etc.: 

Keywords: microscopic nuclear cluster model, generator coordinate method 
Nature of physical problem 

This computer code contructs matrix elements to be used in microscopic studies of light multicluster 
system. 

Method of solution 

We analytically evaluate the determinants and inverses appearing in the matrix elements of Slater deter- 
minants of nonorthogonal single particle states, using the Mathematica symbolic manipulation language. 
A fortran code is used to substitute the space part of the single particle matrix elements. 

Restrictions on the complexity of the problem 
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Varies with available RAM and type of computer. Use of the available RAM depends on the number of 
clusters and nucleons. 

Typical running time 

For the example given in the test run input, approximately 300 s. 
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LONG WRITE-UP 



1 Introduction 

The microscopic nuclear cluster model provides an unquestionably successful description of light nuclear 
systems jlj , Q . The main difficulty in performing microscopic cluster-model calculation is the evaluation 
of the microscopic Hamiltonian between Slater determinants of nonorthogonal single particle-states. This 
calculation becomes especially tedious for multicluster systems. The need for analytical form of these 
matrix elements calls for application of symbolic manipulation languages. 

In this paper the single particle states of the Slater determinants are displaced gaussian functions, 
that is, we determine the matrix elements of the "generator coordinate model" (GCM) M. To achieve 
this, we combine the advantage of computer languages the Fortran and the Mathematica |J: With the 
help of Mathematica we can derive analytical expressions for the matrix elements, with Fortran we will 
sort out and specialize the formulae. 

Our program is applicable for a general N-cluster system of Os clusters of common oscillator width. 
The generalization of the program for nonequal width parameters or for higher harmonic oscillator orbits 
is straightforward. We hope that our program brings the technique of the nuclear cluster model in the 
reach of nuclear physicist. 

2 Formulation of the problem 

Let us consider an iV-cluster system, where each cluster consists of Ai (i = 1,...,N, Yli=i-^i = A) 
nucleons, and the fcth nucleon occupies the single particle state ipj(xk) (see Figure 1). The notation Xk 
stands for the space-, spin- and isospin-coordinate of the fcth nucleon: Xk — Uk, Tfe). We assume that 
the nucleons of the ith cluster occupy the single particle states of the form 

ipj(x k ) = (pi(r k )xa 3 (k)x Tj (k), (j = Ai-i, A,), (1) 
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where ^(r^), Xo-j(k) and Xt<0) are the space, spin and isospin parts of the wavefunction, respectively. 
(Note the cluster label i on ^(r,).) We restrict ourselves to the case, where the space and the isospin 
parts of the single-particle function of the nucleons of the A particle system are fixed, but the spins may 
have different orientation. A given set of spin quantum numbers is labelled by a. 

The spin and the isospin functions of different states are orthogonal 

(Xoklx*,) =<Wn (2) 

(Xr k \Xn) =<Wi, (3) 
the space-part of single-particle states overlap, defining an overlap matrix 

Ki,j) = {<Pi\<Pi)- (4) 

The wave function of the A-particle system is given by linear combination 

* = ^a Q * Q (5) 

a 

of the Slater determinants 

= J=det{^)} (6) 

of the different sets of single particle states. 

The matrix elements of Slater determinants can be calculated using well-known rules |4| , Q . In the 
following we give a concise summary of these formulae to establish the formalism of our program. 

The overlap of two Slater determinants can be written as 

(**|tt*> = det{(^|Vf>}. (7) 
By using eqs. (J2j> , (||) and (|j) the overlap of the single-particle states becomes 

B^ = m\^)} = KhjKa !7? Sr kn , (k = A i - 1 ,...,A i , I A, , (i,j = l,...,N) (8) 

thus to determine the overlap we must calculate the determinant of the matrix 

B af3 = {B«f} U.j ! 1:. (9) 



The matrix elements of the one-body operator are given by 

A A A 

Y.moPl*' 3 ) = ££(CI0 (1) |^)A<f, (io) 



where A"- is the cofactor of B a/3 (i.e. the subdeterminant, obtained by crossing out the ith row and the 
jth column of B a P and multiplying it with the phase (— 

The matrix elements of the two-body operator can be determined from the formula 

A A 

= X](C^|0 (2) I^V'f)a(-l) i+i+fe+ 'detD^, (11) 

i<j i<3 
k<l 

where 

\4>klpl)a = \^klpl) ~ \lpl1pk), (12) 

and where D^ kl is an (A — 2) x (A — 2) matrix built up of elements of B al3 by leaving out its rows i and 
k and its columns j and I. 

In this paper we limit ourselves to the case of spin-isospin independent one-body operators, and for 
the two-body operator of the form 

Oif = V{i, j){w + mP tj r + bPif - hP ZJ T ) (13) 

where are the space- , spin- , and isospin-exchange operators. These restrictions are not too serious, 
one can easily generalize the program to allow for other type of operators. 



3 Mathematica part 

The analytical calculations of the determinants and sums of eqs. (Q), ([h]) and ( pi] ) can be easily performed 
by using the algebraic computer language Mathematica. 

By the aid of Mathematica we first evaluate the single particle overlaps (see eqs. (||) and (H) ) and 
then build up the matrices B, A u and D^ kl and carry out the operations in eqs. (R), ( flO| ) and ( p"l| ) . 
The matrix elements ( "kernels" ) have the form 

- overlap : 

<^i^>=£< n&(^j fe s (u) 

n=l Z„=l 

- one-body : 

A n a b m n 

Y,{* a \oV\*P) = j2c£t(p n>qn ) n &(*, B) i,j*s (15) 

i=l n=l /„ = 1 

where t(p,q) = \<p q ), 



- two-body : 

. tb 

A v tb n„ m„ 

^{y a \0lf\*P)=^(wW u + mM v + bB u + hH v )v( Pv ,q u] r„,s„)J2c% [] 6(< /n> ji n ) fcl -, (16) 

where (<p p „(Pq v \0( 2 > |<£r„¥>s„) = v(p v ,q v ;r u ,s„)(wWv + mM v + bB v + hH v ). All the coefficients, indices, 
limits of summations and products are to be determined by Mathcmatica. A detailed example of these 
formulae can be found in section 3.3. 



3.1 Input data 

The input data are to be read from the input file "nuc.dat". Detailed input data specification: 

1) NumClus: Number of clusters. 

2) NumPar: Number of particles in the clusters. (Numpar(i), i=l,...,Numclus) 

3) isospin: Charge number of the particles, 1 stands for protons and stands for neutrons, (isospin(i), 
i=l,...,A). 

4) SpinConf: Number of spin configurations. 

5) spin: Spin of the particles, 1 stands for spin up, stands for spin down, (spin(ij), i=l,. .., A, j=l,..., SpinConf). 

6) a a : Clebsch-Gordan coefficients of the spin configurations, (a=l,...,SpinConf). 

3.2 Output data 

To calculate the overlap, the one-body and the two-body kernels, one should use the programs "ok.m" , 
"obk.m" and "tbi.m" , respectively. The output is written into the files "oJc.out", "obk.out" and "tbk.out" , 



accordingly. As output we get the coefficients and indices of eqs. (|14|), (JL5|) and (16) ordered according 



to Table 1. In the output files the results for different spin configurations follow in consecutive order. 



3.3 Test run 



As a test run, we consider the example of 6 He—a + n + n. This is a three-cluster system, the two outer 
neutrons may have two different spin configurations with antiparallel spins (see the test run input). The 
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combination coefficients are a\ = —a 2 =< \\\ — ^|00 >= 1/V2- The package must be loaded while 
running the Mathematica program. The program should work with any Mathcmatica implementation 
on any machine without modifications. The output of "ok.m" (ok. out) corresponds to the overlap 

(tf |#) = 6(1, 1) 2 6(1, 2)6(1, 3)6(2, 1)6(3, 1) - 6(1, 1) 3 6(1, 3)6(2, 2)6(3, 1) - 6(1, 1) 3 6(1, 2)6(2, 1)6(3, 3) 

+ 6(1, 1) 4 6(2, 2)6(3, 3) - 6(1, 1) 2 6(1, 2)6(1, 3)6(2, 1)6(3, 1) + 6(1, 1) 3 6(1, 2)6(2, 3)6(3, 1) 

+ 6(1, 1) 3 6(1, 3)6(2, 1)6(3, 2) - 6(1, 1) 4 6(2, 3)6(3, 2) - 6(1, 1) 2 6(1, 2)6(1, 3)6(2, 1)6(3, 1) 

+ 6(1, 1) 3 6(1, 2)6(2, 3)6(3, 1) + 6(1, 1) 3 6(1, 3)6(2, 1)6(3, 2) - 6(1, 1) 4 6(2, 3)6(3, 2) 

+ 6(1, 1) 2 6(1, 2)6(1, 3)6(2, 1)6(3, 1) - 6(1, 1) 3 6(1, 3)6(2, 2)6(3, 1) - 6(1, 1) 3 6(1, 2)6(2, 1)6(3, 3) 

+ 6(1,1) 4 6(2,2)6(3,3). (17) 

There are seven different terms, as the 1st, 5th, 9th, 13th, the 2nd, 24th, the 3rd, 15th, the 4th, 16th, 
the 6th, 10th, the 7th, 11th and the 8th, 12th terms of 16 elements of the sum are equal. These will be 
collected by the Fortran part. 

The output of the programs "obk.m" and "tbk.m" can be understood similarly: 

(*|T|*> - -i(3, 3)6(1, 1) 3 6(1, 2)6(2, 1) + t(3, 1)6(1, 1) 2 6(1, 2)6(1, 3)6(3, 1) + (18) 

and 

= ( Wi + Mi +H 1 +B 1 H1,1 > 1,1)(6(1 > 2)6(1, 3)6(2, 1)6(3,1) (19) 

- 6(1, 1)6(1, 3)6(2, 2)6(3, 1) - 6(1, 1)6(1, 2)6(2, 1)6(3, 3) + 6(1, 1) 2 6(2, 2)6(3, 3)) + 
( W 2 + M 2 + H 2 + B 2 )v(l, 1, 1, 1)(6(1, 2)6(1, 3)6(2, 1)6(3, 1) 

- 6(1, 1)6(1, 3)6(2, 2)6(3, 1) - 6(1, 1)6(1, 2)6(2, 1)6(3, 3) + 6(1, 1) 2 6(2, 2)6(3, 3)) + .... 

As the output of "tbk.m" is too long (560 terms), we present only the first 546 lines of "tbk.out" . 

The "-1" at the bottom of the Mathematica's output serves as an end of file indicator for the Fortran 

part. 



4 Fortran part 



Due to the spin-isospin symmetry, a great number of terms of the sums (|l4|), ( [l5|) and ( |l6|) will be exactly 
equal in the Mathematica's output. By collecting these, the number of terms can be reduced by about an 
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order of magnitude. E.g., in the case of s He=a+n+n+n+n five cluster system, the potential energy kernel 
contains about 50000 terms (see eq.(16)),and only 3528 of them are different. As the Mathematica stores 
the symbols and operations, this collection by Mathematica would require enormously large memory 
and might be very slow. As in the Mathematica's output the analytical formulae (|l4|), ( [l5|) and ( |l6| ) 
are represented by integer numbers, this simplification can efficiently be achieved by using the Fortran 
language. 

After collecting the identical terms, the forms of the expressions (|TJ), ( |f5| ) and (|l^) remain the 
same, except for the coefficient c° , c° b and , where the multiplicity of each different terms appears. 

In the Fortran part we can also substitute the concrete form of the spatial part of the overlap, 
one-body and two-body matrix elements of the single-particle functions. In this paper we limit ourselves 
to the case of displaced gaussian single-particle functions of equal size parameters: 

^)=(v) e ~' y(r ~ Sl)2 ' ( 2 °) 

where is the "generator coordinate vector" pointing to the centre-of-mass of the zth cluster (see Figure 
1). For the case of nonequal size parameters or for other type of single particle functions this part of the 
program can easily be generalized. As the generator coordinates in the bra may differ from those in the 
ket the vectors in the bra and ket position will be distinguished by using a prime ('). 

The two-body potential is represented by gaussian form factor 

V(i,j) = e-*«-'tf ) (21) 

and as an example, the one-body operator is chosen to be the kinetic energy operator T. The single 
particle matrix elements are [gj 

b(i,j) = e<«-^', (22) 
h 2 

t(p, q) = b(p, q) — v{Z - v(s p - s',) 2 ), (23) 
2m 

v(p,q;r,t) = 6 3 / 2 b(p,r)b{q,t)e- 1 ^ 1 ^ s ''- s '' +s ' r - s ' t)2 , (24) 

where = v j(y + /i). 

By substituting the single particle matrix elements (p2|), ( pj ) and (|24|) into eqs. (|l4|), ( |l5| ) and ( |l6| ) 
the matrix elements of the wave function ^ take the form: 
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- overlap: 

(*|*)=£<7°exp(--^U° s ), (25) 
i=i z 

- one-body: 

<*I T I*> = £ £ ^ " ^^f s)exp(-^ S t^6 fl)) (26) 

- two-body: 

N t b -, 1 — 

(flog) |*) = 3 / 2 £>Cf + mCt 1 + bCf + W7?)exp(--i/*t^ s - v (a?s) 2 ), (27) 

i=l 

where A°, Af, Af and £?° b are 2N x 27V matrices, af is an 2AT dimensional vector, and 

= (S1,S2, ...,Sjv,s'i,s'2, ...,s'jv) • (28) 

After suitable angular momentum projection /citeBrink, these matrix elements are building blocks of the 
cluster model calculations. 

4.1 Input data 

The Fortran part needs the input data stored in the file "nuc.dat", and the Mathematica's output 
" ok. out" , "obk.out" and "tbk.out" . Free format is used. 

4.2 Output data 

The Fortran programs "oic.f', "obk.F and "tbk.f handle the overlap, the kinetic energy and the potential 
energy kernels, respectively. The output is written in the files "oi.inp", "obk.inp" and "tbk.inp" ordered 
as shown on Table 2. 

4.3 Test run 

After collecting the identical terms, the Fortran part substitutes the single particle matrix elements, thus 
e. g. the overlap looks like: 

(tf = 2cxp{-^(2(s 1 - s\) 2 + (si - s' 2 ) 2 ( Sl - s' 3 ) 2 + (s 2 - s'O 2 + (s 3 - s'i) 2 )} 
- exp{-^(3(s 1 - s'x) 2 + (si - s' 3 ) 2 + (s 2 - s' 2 ) 2 + (s 3 - s'i) 2 )} 
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- exp{--(-3( Sl - s\) 2 + (si - s' 2 ) 2 + (s 2 - s'x) 2 + (s 3 - s' 3 ) 2 )} 
+ exp{-^( + 4( Sl - s'x) 2 + (sa - s' 2 ) 2 + (s 3 - s' 3 ) 2 )} 

exp{-~(3(si - s'x) 2 + (81 - s' 2 ) 2 + (s 2 - s'3) 2 + (s 3 - s'O 2 )} 

- exp{-~(3(si - s'x) 2 + (si - s'3) 2 + (s 2 - s'O 2 + (s 3 - s' 2 ) 2 )} 
+ exp{-^(( Sl - s\) 2 + (s 2 - s'3) 2 + (s 3 - s' 2 ) 2 )} 



From this form the matrices A° in eq. (25) are determined and written into " ok. out" (see the test run 
output). The output files of "obJc.out" and "obk.out" are too long, and only the first 53 and 57 lines are 
listed. 
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Tabic 1 



overlap one — body two — body 

a,p,n a,f3,n ob a, (5 

mi pi,qi W 1 ,M 1 ,B 1 ,H 1 

ill, ill, fcii mi pi,q!,n,si 

■ iii, Aii n i 

imi , imi , A mi . 777.1 

c l imij i imij , k mil *li, ilu All 

m 2 cf : 

il2,il2,Al 2 P2,<l2 imn jmi> k mi 

: m 2 cf 

: iii, iii, Aii 

im 2 , im 2 , A m2 
„t6 



iii, iii, Aii 



im ni ? im ni , A TOni 



„tfc 
Wii 



a,/? 

P2, 92,7-2,52 

n 2 ?77i 
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Table 2 



overlap one — body two — body 

N N ob N tb 

Co /~iob r<W riM /~<B /~iH 

(Af) n ... (A?) 1JV «) n ... «) 1JV (Af) n ... (A? 



- (Bf) NN 



r<w r<M r<B r<H 

°2 ' 2 i°2 i °2 



liV 



C°2 «) n ••• (B?) 1N «) x (af) N 
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